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Abstract 

We propose a system of nonlinear integral equations (NLIE) which 
gives the free energy of the Uq{sl{r + l|s + 1)) Perk-Schultz model. In 
contrast with traditional thermodynamic Bethe ansatz equations, our 
NLIE contain only r + s + 1 unknown functions. In deriving the NLIE, 
the quantum (supersymmetric) Jacobi-Trudi and Giambelli formula 
and a duality for an auxiliary function play important roles. By using 
our NLIE, we also calculate the high temperature expansion of the free 
energy. General formulae of the coefficients with respect to arbitrarily 
rank r + s + 1, chemical potentials {^a} and q have been written down 
in terms of characters up to the order of 5. In particular for specific 
values of the parameters, we have calculated the high temperature 
expansion of the specific heat up to the order of 40. 

MSC: 82B23; 45G15; 82B20; 17B80 

PACS2003: 02.30.Rz; 02.30.Ik; 02.20.Uw; 05.70.-a 

Key words: nonlinear integral equation; Perk-Schultz model; quantum Jacobi- 
Trudi and Giambelli formula; quantum transfer matrix; thermodynamic Bethe 
ansatz; T-system 
OIQP-05-12 

1 Introduction 

In statistical physics, to calculate the free energy of solvable lattice models 
for finite temperature is one of the important problems. For this purpose, 
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thermodynamic Bethe ansatz (TBA) equations have been often used p. 
In general, the TBA equations are an infinite number of coupled nonlinear 
integral equations (NLIE) with an infinite number unknown functions. Then 
it is desirable to reduce TBA equations to a finite number of coupled NLIE 
with a finite number of unknown functions. 

Destri, de Vega |2] and Kliimper [21 E] proposed NLIE with two (or one 
^) unknown functions for the XXZ (or XYZ) model. To generalize their 
NLIE to models whose underlying algebras have arbitrary rank seems to be 
a difficult problem as we need considerable trial and errors to find auxiliary 
functions which are needed to derive the NLIE. Then there are NLIE of 
abovementioned type for models whose underlying algebras have at most 
rank 3 (for example, 00111). 

Several years ago, Takahashi discovered jH] an another NLIE for the XXZ 
model in simplifying TBA equations. Later, the same NLIE was rederived 
jH] from fusion relations (T-system) among quantum transfer matrices 
(QTM) jTTj. In addition, it was also rederived ^12j for the XXX model from 
a fugacity expansion formula. 

In view of these situations, we have derived NLIE of Takahashi type for 
the osp{l\2s) model the sl{r + 1) model [T^, the higher spin Heisen- 
berg model [121, the Uq{sl{r + 1)) Perk-Schultz model [15]. In these cases, 
the number of unknown functions and NLIE coincide with the rank of the 
underlying algebras. In this paper, we will further derive NLIE with a finite 
number of unknown functions for the Uq{sl{r + l|s + 1)) Perk-Schultz model 
pri llBj. which is a multicomponent generalization of the 6- vertex model and 
one of the fundamental solvable lattice models in statistical mechanics. For 
example, a special case of this model is related to the supersymmetric t — J 
model, which is important in strongly correlated electron systems. 

In section 2, we introduce the Uq{sl{r + l\s + 1)) Perk-Schultz model, 
and define the QTM for it. As a summation over tableaux labeled by a x m 
Young (super) diagram, we introduce an auxiliary function ()2.13|) 1201 
which includes an eigenvalue formula ()2.10p of the QTM as a special case. We 
also introduce a system of functional relations (T-system) which is satisfied 
by this auxiliary function. 

In section 3, we derive two kind of NLIE which contain only r + s + 1 
unknown functions. The first ones ()3.9p . ()3.10|) reduce to the NLIE for the 
Uq{sl{r + 1)) Perk-Schultz model in ||T2] if s = —1. However our new NLIE 
are not straightforward generalization of the ones in our previous paper JH] • 
In fact for r, s > straightforward generation of our previous NLIE 

becomes a system of an infinite number of coupled NLIE which contains an 

^if an integral contour with a closed loop is adopted 
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infinite number of unknown functions (see ()3.7p ). To overcome tliis difficulty, 
we will use the quantum (supersymmetric) Jacobi-Trudi and Giambelli for- 
mula ()2.25|1 and a duality ()2.23|1 for the auxiliary function, from which a 
closed set of NLIE can be derived. We will also propose another NLIE ()3.2H) 
and ()3.22|1 in the latter part of the section 3, which have never been con- 
sidered before even for the Uq{sl{2)) case. In deriving the NLIE, we assume 
that q is generic. However we expect that our results can be also analytically 
continued to the case where g is a root of unity. 

In section 4, we calculate the high temperature expansion of the free 
energy based on our NLIE. In particular, we can derive coefficients fl4.2|) - ()4.6|) 
up to the order of 5 for the arbitrary rank r + s + 1. The point is that if we fix 
the degree of the high temperature expansion, we can write down a general 
formula of the coefficients. On the other hand, if we specialize parameters, we 
can derive the coefficients for much higher orders. For example for (r, s) = 
(2, — 1), (— 1, 2), g = 1, = case, coefficients of the high temperature 
expansion of the specific heat up to the order of 40 are presented in appendix. 
It will be difficult to derive the coefficients of such a high order by other 
method. 

Section 5 is devoted to concluding remarks. 

2 The Perk-Schultz model and the quantum 
transfer matrix method 

In this section, we will introduce the Uq{sl{r + l\s + 1)) Perk-Schultz model ^ 
da HE] and the quantum transfer matrix (QTM) method [11112212311231113 

for it. The QTM method was applied to the Perk-Schultz model in ref. |Sj 
(see also, ref. |2Sl EEl El ) • 

Let us introduce three sets B = {1,2, ...,r + s + 2} = U B_, where 
5+ n 5_ = 0, = r + 1 and = s + 1 (r, s G Z>_i). We define a 
grading parameter p{a) (a e B) such that p{a) = for a G -B+ and p{a) = 1 
for a G B_. The i?-matrix of the Uq{sl{r -|- l|s + 1)) Perk-Schultz model dZj 
is given as 

R{v) = Yl Kli'ME^''"' ® E"^'''^ (2-1) 

ai,a2,f)i,b2G-B 

where i?"''' isar + s + 2byr + s + 2 matrix whose element is given as 

'^Uq{sl{r + l|s + 1)) is a quantum afRne superalgebra, which characterizes the i?-matrix 
of this model. See for example, |^. We assume 77 S K (g = e"^). A rational limit {q 1) 
of the Perk-Schultz model is the Uimin-Sutherland model [271 128j . 
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{E""'')ij = 5ai5bj] Klili'") is defined as 

K:>) = [{-iy^% + 1]„ (2.2) 
<(^) = (-ir^"^^^'^H. i^^b), (2.3) 

where G C is the spectral parameter; a,b E B; [v]q = (g^ — q~^)/{q — q^^)] 
q = e''. Note that this i?-matrix reduces to the one for the well known 
6- vertex model if (r, s) = (1, —1). 

Let L be a positive integer (the number of lattice sites). The row-to- row 
transfer matrix on (^([^r+s+2^(»L jg defined as ^ 

tiv) = tr-oiRoLiv) ■ ■■Ro2iv)Roi{v)). (2.5) 

The main part of the Hamiltonian is proportional to the logarithmic deriva- 
tive of the row-to-row transfer matrix ()2.5j) : 

^ ^ j=l ^ a&B 

J2 (sign(a - b) sinhr^E;''^^;;^ + (-1)p('^M>^) Ef'' Ej^,) \, (2.6) 

a,beB 

where we adopt the periodic boundary condition E'^^^ = E^'^ . Without 
breaking the integrability, we can also add the chemical potential term 

L 

Hck = -Y.Y.f^^r (2-7) 

j=l a£B 

to Hbody Then the total Hamiltonian is H = Hbody + Hch- 

To treat the model at finite temperature T, we introduce the so-called 
quantum transfer matrix (QTM)^^: 

tQTuiv) = Yl ^QTM(^){S:.f4^f"^^2^" ■ ■ ■ ^S'""' (2-8) 

JV 

■^The lower index i,j of Rij{v) is used as follows: for example, E"^''^ is defined on 
((|-r+s+2)®(L+i). ^a,b ^ jm ^^a.b ^j(g,{L~k) ^ ^^^^^ j ^^^^2 by r+s+2 identity matrix; 

fc = 0,l,...,L. Theni?„(^;) is defined as % ^ = Ea,,a.,6„6. K^;^i(^')£^r''''£^'"'''- The 
trace tro is taken over the auxiliary space indexed by 0. 
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where N G 2Z>i is the Trotter number; Un+i = i^i', ^kiCik^Pk E B; u = 
-^H^; Klili^) = Rbl'Z^'"^ is '90° rotation' of R{v). We can express 
jTTj the free energy per site in terms of only the largest eigenvalue Ai of the 
QTM at = 0: 



/ = -T hm logAi, (2.9) 

where the Boltzmann constant is set to 1. 

Due to the Yang-Baxter equation, the QTM ()2.8|) forms commuting fam- 
ily for any v. Thus it can be diagonalized by the Bethe ansatz. The eigenvalue 
formula ^ of the QTM will be (cf. 

T«(^;)=^^(a;^;), (2.10) 

where 



z{a;v) = iJa{v)^a 

ft_.(„- -sa^-"- )Q„(.-!j;;=^) 

Ma 

Qa{v) = llsmr]iv-vl;^), (2.11) 

k=l 

N 

'sin r](v ± iu) 



sinh?7 



where Ma G Z>o; Qo{v) = Qr+s+2{v) = 1. ^a£{ — l,l}isa parameter which 
depends on the grading parameter {p{b)}b(zB- {^1^^} is a root of the Bethe 

'^To be precise, this formula is a conjecture for general parameters r, s,q, Ha, N . In 
0, the algebraic Bethe ansatz for a one particle state was executed for the QTM of the 
Uq{sl{r + l\s + 1)) Perk-Schultz model. As for the Uq{sl{2)) case, a proof of this formula 
by the algebraic Bethe ansatz is similar to the row-to-row transfer matrix case (cf. j29|). 
This formula has a quite natural form (dressed vacuum form) from a point of view of the 
analytic Bethe ansatz PUI |^ ■ An eigenvalue formula of the row to row transfer matrix 
H2.5() was derived in |32l [TH] . It has essentially same form as H2.10|l except for a part which 
is related to the vacuum eigenvalue. There is also support by numerical calculations for 
small r, s. 
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ansatz equation (BAE) 

V'aH-i(4'^^ + tE;=i(-iM^-)) 

for A; G {1,2, . . . ,M„} and a G {1, 2, . . . , r + s + 1}, 

where Sa = G !}• From now on, we assume the relation p(l) = p(r+ 
s + 2) on the grading parameter. In this case, the eigenvalue formula ()2.1()|1 
of the QTM has good analyticity to derive the NLIE. We expect that this 
assumption does not spoil generality as the free energy will be independent 
of the order of the grading parameters. 

Let us define an auxiliary function [121120] (see also [H^): 

a m 

T^\^)= 5^ nn^(^^>'^-^(«-^-2j + 2fc)), (2.13) 

where m,a E Z>i, and the summation is taken over dj^k G -B (1 < 2 < ■ • • < 
r + s + 2) such that 

dj,k < dj+i^k and dj^k < dj^k+i (2.14) 
dj,k < dj^k+i if dj^k G or dj^k+i e (2.15) 
dj^k < dj+i^k if dj^k e B+ or dj+i^k e B+. (2.16) 

This function contains tI'^\v) ()2.1()p as a special case (a, m) = (1, 1). ()2.13|1 
can be interpreted as a summation over a Young (super) tableaux labeled 
by a X m Young (super) diagram. It is related to a system of eigenvalue 
formulae of the QTM for fusion models [31]. Note that the condition ()2.15|) 
is void if s = —1, then ()2.13|) reduces to the Bazhanov-Reshetikhin formula 

insj. 

For a,m E Z>i, we will normalize ()2.13j) as Tm\v) = Tm\v) / N'm\v) , 
where 

X n n '^-(^ ^ — ^ — 0- (2.17) 

j=i k=i 
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Here we introduce a parameter ^ G { — 1, 1}. Tm\v) has no pole on v due to 
the BAE dsn. In contrast, f^\v) has poles at v = ±{^^i + iu) + ^ 
(n e Z) for (a, m) G Z>]^x {1, 2, . . . , s + 1} U {1, 2, . . . , r + 1} x Z>i. 

One can show that Tm\v) satisfies the so called T-system for Ug{sl{r + 
l\s + 1)) jini ED] (see also [SE! for a derivation of TBA equations from the 
T-system). For m, a G Z>i, 

TL^Kv - '^)f!^\v +'-)= fi:Uv)Tl%M + f^^-'Kv)f^^^'\v) (2.18) 
for a G {1, 2, . . . , r} or m G {1, 2, . . . , s} or (a, m) = (r + 1, s + 1), 

Ti:^'\v-l)f^Z^'\v + ^)=f'^ll\v)f^^^^^ for mGZ>,+2, (2.19) 
Ti:\(^;-^)ri:U^^ + ^) = Ti:T^^(^)Ti:r^(^^) for aGZ>,.+2, (2.20) 



where 



(t; + ?z)0+(t;-^z) 



There is a duality relation for the auxiliary function. 

fi;7H^)=C-'Ti;r^(^) for aeZ>i, (2.23) 

Ma 

where C = + ^"'X ■ (IT^ (resp. ^TT^) becomes 1 if { = 1 (resp. 

^ = —1). Note that there is no upper bound for the index a of Tm\v) 
for m G {1, 2, . . . , s + 1} if s G Z>o. For s = —1, this T-system reduces 
the one for Uq{sl{r + 1)) jSEl (see also [Ej). In this case, ()2.23|1 reduces to 

fi:t'V) = C"' = e '°"''^'^'^T^'"^"'"''' if e = 1 (see eq. (2.21) in [12]). From 
the relations (jTTT^ . (jT^ . (jT^ and pTR|l for (a, m) = (r + 1, s + 1), one 
can derive the following relation for a G Z>2: 

Ti;rH^) = c-^Ti;r^(^) 

n;=2(cTr'^(^ + + t^%{v + ^^)) ■ 

(2.24) 
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Tm\v) can also be written in terms of a determinant (the quantum (super- 
symmetric) Jacobi-Trudi and Giambelli formula fH^ ED] (for s = —1 case, 
ISni; for UgiBi^'^) case, |SII)) 

- ^^"^(^)JJ,.(^^V.(-^^^^))^ (2.26) 

where fi''\v) = for a < and fjn\v) = for m < 0. W^\v) and zt\v) 
are normalization functions: 

W^\v) = — , (2.27) 



2 

where 11^=1 (' ■■) = !■ Substituting p.25|l into p.23|) . we obtain an equation 

l<j,fc<a+s \ \ Z 

= C-'wilT\v) det ('fr'^+^-'=) ft 

l<j,fc<s+l Y y 

for a G Z>i. (2.29) 
Expanding partially ()2.29|) on both side, we obtain 



J + A; - s - 2 

f z 

2 



{-ir^^A,{v) + (-i)^e-^S^yS^4(t;) 

for a e Z>2, (2.30) 
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where 



M^) - , Jt,. (a. - i±i^|^.) ) (2.31) 

fi.tM^Ti'-+'*'-'''(v) for {j,k)ji(a + s,l), /„+.,,(«)= 0, 

9MM=f<'+"+'-''">(i,) for + g.+uW-O, 

A,{v)= det |^^_ ^' + ^-'^-\ ^y (2.34) 



It turns out that T]^"^'"^*'*(t') is written in terms of {tI^^\v)} where max(0, r — 
s + 2-a) <d< a + r + s-l. Then fi^"^(t;) for a G 'L>r+s+2 can be expressed 
in terms of {t['^\v)} where Q<d<r + s + l. Similarly, we can derive the 
following relation from ()2.23|) and ()2.26|) . 

for a e Z>2, (2.35) 

where 

T/N , (~r ( a + r + 1 — ?' — A;\\ 

^=<"' = .^..r K ^"^ V ) 

hj^kiv) = T;+i+j._^(i;) for (j, A;) 7^ (a + r, 1), /ia+r,i(t') = 0, 
M") = for (j, + 1,1), 6,+„(«) = 0, 

Let us consider the limit 

a m 
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where the summation is taken over {dj^k} {dj,k £ B) which obey the rules 
(ini-dnni)- For example, for Uq(sl{2\l)) (5+ = {1,3}, B_ = {2}) case, we 
have, 

Q? = eie- +66^+^6-, (2.41) 

= ^2"^e^^^Qf^ for aGZ>2. (2.42) 
We can also rewrite ()2.42p as 

= ^ = — ^^T^ . (2.43) 

This quantity ()2.4Up corresponds to the character of a-th anti- (super) symmetric 
and m-th (super) symmetric tensor representation. We will use Q^^^ and Qm'' 
later. 

Qm^ also satisfies the so called Q-system, which is the T-system p.l8|l - 
()2.23|1 without the spectral parameter v. for m,a & Z>i, we have 

QL"^' = QtliQtli + Qt'^Qt^'^ (2.44) 

for a G {1, 2, . . . , r} or m G {1, 2, . . . , s} 

or (a, m) = (r + 1, s + 1), 

Q!;^'^' = for mGZ>,+2, (2.45) 

Qii' = QiVi^Qiti^ for aGZ>,+2, (2.46) 



where 



Qo"^ = Q^^ = l for a,mGZ>i, 

= C-^Q^+r^ for aGZ>i. (2.47) 



The Q-system was introduced [o8j 122] as functional relations among char- 
acters of finite dimensional representations of Yangians (or quantum affine 
algebras) associated with simple Lie algebras. The above system of equations 
is a superalgebra version of them. 

In closing this section, let us comment on the analyticity of the auxiliary 
function ()2.13|) . As mentioned before, the free energy ()2.9|) is given only by 
the largest eigenvalue of the QTM ()2.8p . Then we are only interested in 
a root of the BAE ()2.12j) which gives the largest eigenvalue of the QTM. 
Judging from numerical calculations |2S1 EEl El CEl , such a root will exist in 
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the sector ^ = Mi = ■ ■ ■ = M^+s+i of the BAE, and it will form 'one-string' 
on the complex plane. For this root, the zeros of the auxiliary function ()2.13|) 
will exist near the lines Imv = i^^^ at least for {^a} = {0} and small u 
(see, figures in [201 CH CEj)- In this sector, we have 

a = 1 for be B, 

eb = l for beB~{r + s + 2}, (2.48) 
C = exp( ). 

From now on, we only consider the largest eigenvalue of the QTM, and 
assume these values ()2.48j) of the parameters. 



3 The nonlinear integral equations 

In this section, we will derive NLIE by using formulae in the previous section. 
We will treat two kind of NLIE paying attention to the value of the parameter 
^ G {—1, 1}. Although the first step of calculations ()3.1|l - ()3.fij) is similar to 
s = — 1 case PIEIIISI; we will present it for reader's convenience. 

Taking note on the limit ()2.40|) and the fact that Tm\v) has poles at 
V = ±i^^i + iu) + ^ {ne Z) for (a, m) G {1, 2, . . . , r + 1} x Z>i U Z>i x 

{1, 2, . . . , s + 1}, we can expand Tm\v) as follows. 




(3.1) 



m,j ^ mj 



where the coefficients A^^^-,A^^^- G C can be expressed as contour integrals: 

a'::, = + =^f» + (3.2) 

Here the contour Cm (resp. ) is a counterclockwise closed loop which 
surrounds v = ^-Y^^i+iu (resp. v = —^^^^i — iu) and does not surround v = 
-^^i-iu-^, ^^i+iu+f, (resp. v = -^^i-iu-f), 
where n e 1^, k e "Z — {0}. Using the T-system ()2.18|) - ()2.20|) . we can rewrite 
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dsn as 

M ^ I dt; [ gi(t;-|)fit(^-|) 
_ / d.; r gi(^ + f)gii(^^ + f) 

- few ^ ~ ^^-^^ 



^ — ^i + iuy , 



where we admit fi^V) = if {h,n) e Z>^+2 x Z>,+2 (cf. IHlllII). Substi- 
tuting ()3.3p into ()3.1|) and taking the summation over j, we obtain 

1 



dy 



X 



(a) ( I a+m~l f • , „-„A'r(") 



^ ^(«) 2-Kiv -y - - iu 



TrZi(y + ^^e^ + m)Tr|i(|/ + ^^e^ + lu) 
^ f^-'\y + ^±f^e^ + iu)fl^^'\y + ^^^^ + lu) 



1 



N_ 
2 



2mv - y + + iu + ^ 



c. 



I T^^\y - ^^ii - lu) 

^ fi^~'\y - ^^e^ - tu)ft^'\y - ^^^ji - lu) 
f^\y - ^^^^^t - tu) 

Here the contours are shifted as follows: the contour Cm (resp. C„ ) is a 
counterclockwise closed loop which surrounds y = (resp. y = 0) and does 
not surround?/ = —{a+'m)C,i—2iu— — , — (resp. y = {a+m)C,i+2iu+ — , —), 
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JV 
2 



where n E Z, k E Z — {0} . We can neglect the terms ( , a+m/, ■ , ,rn I in 

\ 2 ^ >7 y 

()H.4j) since the poles at ?/ = in the two brackets {■ ■ ■ } are canceled by the 
zeros from these terms. By using the following relation 

m _ 

lim y -—- = (3.5) 
m^oo ^-^ V — — tan™ 

n=-m V ' 

we can take the summation over n E Z. 
+ 



c(^) 2m tan r]{v - y - - zn) 

fj^Uy + ^^e^ + ^u)f;^Uy + ^^^^ + ^u] 



^ Tt-'\y+ ^^e^ + zn)Tir+'^(y + ^^^^ + m; 



+ 



dy 7] 



gj^' 27ri tan f]{v - y + + to) 



- - m)Tirji(y - - m) 

I T4")(i/-^^e^-2w) 

f{^~l)^y _ a+r^^, _ ,^)fir+^^(l/ - ^^jl - lu) \ 

+ fir)(i/-^e*-m) J' 

for (a, m) G {1, 2, . . . , r + 1} X Z>i U Z>i X {1, 2, . . . , s + 1}. 

In the next subsection, we will consider specializations of this system of NLIE 

(EHD. 



3.1 The nonlinear integral equations for ^ = 1 

Let us consider the NLIE ()3.6j) for ^ = 1 and m = 1. Taking note on the 
fact f^"'\v) = 1 (cf.(I22H)), we can drop the first terms in the two brackets 
{■ • • } in ()3.6|) since they have no poles at ?/ = 0. Then the NLIE ()3.6|) reduce 
to the following NLIE on T^'^\v) = limAr__^oo ^i"'' ("w) after the Trotter limit 
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N ^ oo with u 



J sinh rj 
riNT ■ 



+ 



dy yrr^'iy + Wiiy + f) 

for a G Z>i, (3.7) 



where the contour c|"^ (resp. C^""*) is a counterclockwise closed loop around 
y = (resp. y = 0) which satisfies the condition y ^ v — + ^ (resp. 

y ^v + ^i + ^) and does not surround zi"^ - + ^' + + T' T 

(resp. + ^2 + ^, (a + l)i + ^, (n G Z, A; G Z - {0}). Here we put 

the zeros of (f ) as (ti) is a known function: 

^°''")->"^°''")--p( r(co":^'''c:l„)) V (-) 



Note that ()3.7|) are an infinite number of couple NLIE if s G Z>o. This 
situation is quite different from the Uq{sl{r + 1)) case [121 Ellin]- However 
these NLIE are not independent, then we will take the first r + s + 1 of them 
(dSHI) for a G {1, 2, . . . r + s + 1}). The NLIE for a = r + s + 1 contains 
Ti^^^'^\v), then we will eliminate this by using the relation ()2.3()|1 . where 
W^\v) = 1 for e = 1. 



T/")(^;)=Q(") + 



dy vTl''-'\y + f)Tl^^'\y + f) 
27n tan 77(1; -y- ^ll)rl''\y + 

ia 
2 



+ 



ci"' 27r« tan ti{v - y + ^11)tI^^ (y - 



for a G {1,2, ...r + s}, (3.9) 

T/'■+^+l)(^;)=gf+^+l) 

dy vTl'-^'Hy + '-^^)Hy + 



+ 



^(.+^+1) 27r2 tanr/(t; -y- l(l^)r/''+^+^)(y + 



2 vy ' 2 

-(r+s)/ 1(!±^±1) ^ T/'., j(r+£+l)^ 



+ / dy r^rr'\y-'-^)Hy- , . .3 

7cr^+i) 27rz tanr7(t; - y + !(ld^)T/'^+^+i)(^ _ l(l+£)) ' ^ ' ^ 

^(.) = hm fr-^)(.) = , (3-11) 
^' ^-°o ' (-l)«A3(t;) + (-1)-C^4(^) ^ ' 
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where 

f,.l.(v) = T!'*'+'-''\v) for U,k)^(s + 2,1), = 0, 

9,.t{v)=T;'"*'-'-\v) tor {j.k)=l{s + l,l), g,+u(«) = 0, 

If s = —1, then Ai{v) = A^^v) = and A2{v) = A^i^v) = 1, and consequently 
reduces to J-'{v) = Ti'^^^\v) = Q^i'^^'^ = C = ^~ 5^-^^^ where the 
determinants should be interpreted as deti<j_jfc<o(- ■ ■ ) = 1, deti<j fc<_i(- ■ ■ ) = 
0. Thus dSH) and (imUl reduce to the NLIE for Ug{sl{r + 1)) in In 
particular for s = {Uq{sl{r + 1|1)) case, we can use p.24|) : 



c'"' 27r2 tanr/(t; -y- i^ll)Tl^\y + 
dy r^rl''-'\y - f)T}^^'\y - f) 



for a G {l,2,...r}, (3.16) 



+ 



dy vrl'\y + ^^)Tr'\y+'^: 



+ 



cr+^' 27r2 tan?7(t; - y - ^)(C + T}^'\y + ^)) 

cr^' 27r^ tanr/(t; - y + ^)(^ + r/^)(^ _ iilfll)) " 

(3.17) 

The free energy per site is given by a solution of these NLIE (j3:9|l - (j3T7j) 

/ = Jcoshr/-Tlogr/^\0). (3.18) 

In these NLIE fl3.9|l - ()3.17|l . the number of unknown functions and equations 
is r + s + 1, which contrasts with TBA equations 02 ESI El EEl ES] • 
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3.2 The nonlinear integral equations for ^ = — 1 

Next, let us consider the NLIE ()3.6|) for ^ = — 1 and a = 1. Taking note 
on the fact Tm\v) = 1 (cf. ()2.22|) ). we can drop the second terms in the two 
brackets {■ ■ ■ } in since they have no poles at y = 0. Then the NLIE 
(jH.fij) reduce to the following NLIE on Tm\v) = limN^ooTm\v) after the 
Trotter hmit ^ oo with u = 



i{m—l) •• 



+ 



2 



2 

for m G Z>i, (3.19) 



where 



and the contour (resp. C ^ ) is a counterclockwise closed loop around 
?/ = (resp. y = 0) which satisfies the condition y ^ v + ^^^^i + ^ (resp. 

y^v- + ^) and does not surround + + ^, (1 + m)z + ^, 

f (resp. - + ^, -(1 + m)2 + ^, ^) (n G Z, A; G Z - {0}). Here 

{zm''} are zeros of Tm\v): Tm\zm) = 0. These are an infinite number of 
coupled NLIE. We can reduce them as ^ = 1 case. By using ()2.35|) in the 
limit oo, we can reduce fj3.19|) as follows, where Zm\v) = 1 for ^ = — 1. 

ra)(.) = Qa) + / yT^l^y - ^f)7^Uy - ^f) 



+ 



2 Vtf 2 ^ 

dy r/ri!j,(y + f)ri^_;,(y + f) 

27rz tanr/(t; - y - '^)Ti'\y + ^) 

for m G {l,2,...r + (3.21) 
M = na) , / d, vT^i^^jy - ^)g(y - ^) 
^"^"^^^ ^^"^"^ %..2.nan,(.-, + i(r^)T(^Ui(l/-^^ 

di/ vT}Hiy + '^^)Qiy + '^) (3^22) 



c<t,, 2vrz tanr/(^; - y - ^^)T«+,(y + ^ 
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Qiv) = lim T;'„,Jv) = , ^ , , , (3.23) 

where 



AM = det ( h,,, { V - + ^ J ^ ) ) (3.24) 



l<j,k<r+2 
-(1) 



= for (j,A;)^(2 + r,l), K+2,i{v)=0, 

Aiv) = ^Jet^^^ (b,, (. - ^ + ) (3.25) 

hkiv) = T}X^_,{v) for (j,A;)^(r + l,l), 6,+i,i(t;) = 0, 

r + 2 — j — 

= 1 JS'k. l'^" - '-^K^' I I . (3.27) 

where Tm\v) = for m < 0. 

In particular for r = {Ug{sl{l\s + 1)) case, we can use ()2.24|1 : 



Mv) = , . det , ^ r«,_, - ^ ) ) , (3.26) 



"* (f) Qm"* 



2 ^ 



c<„^) 27rz tanr/(t; - y + l(!fi))ri^)(y - 



2 y^'" 2 
-d) , ^m^-^(l) 



dy r^rj,!^i(y + f)rj,^i(y + f^) 



c'l^ 27rz tanr7(t; - y - l(!Ifll)T^i)(y + li^^; 
for m G {l,2,...s}, (3.28) 

dy r,r«(,-i(fi))r«(,-ii^) 



+ 



+ 



c« 2vrz tanr^(t; - y + ^)(C-i + Ti^)(y - ^)) 
r,T«(y + ^)ri_;}(y + ^) 



/c«, 2vrz tan viv - y - ^)(C-i + r;^)(2/ + ^)) ' 

(3.29) 

The free energy per site is given by a solution of these NLIE ()3.2H) - ()3.29|1 
/ = -Jcoshr7-Tlogr/^^(0). (3.30) 

In some sense, these NLIE are 'dual' to the ones in the previous section. The 
NLIE ()3.21|) - ()3.29p have only r + s + 1 unknown functions. These NLIE have 
never been considered before even for Uq{sl{2)) case. 
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4 High temperature expansions 

In this section, we will calculate the high temperature expansion of the free 
energy from our new NLIE. For large T/| J|, we assume the following expan- 
sion : 

r,(^)(^)=exp 5^en^)(^r+o((^) 

\n=0 

= Q? 1 1 + i>'t\vd + Lptv) + dr + (4.1) 




where b^o\v) = logQi"^- Here we do not expand {Qf'''}b>i with respect to ^. 

Thus the coefficients {bn^ (v)} themselves depend on ^. In this sense, our high 
temperature expansion formula is different from ordinary one. Substituting 
this ()4.H) into some of the NLIE ()3.7p - ()3.17p . we can calculate the coefficients 
{bn\v)} up to the order of n = deg. Note that we only need {bn\o)} to 
calculate the free energy ()3.18|) . Taking note on this fact, firstly we use 
^ a subset (NLIE for a G {1, 2, . . . , de^}) of the non-reduced NLIE (jSHI) 
rather than the reduced NLIE We have observed that b^n\o) 

can be expressed in terms of ^ Q^i\ Qf\ . . . , Q^^~^^\ We have calculated 
the coefficients by using Mathematica. As examples, we shall enumerate the 

^As for numerical calculations of the free energy, we expect that the reduced NLIE 
(|3.9() - (|3.17() are easier to use than the non-reduced NLIE (|3.7|) . 

^For s — —1 case, they are Q'"i\q'"i \ • ■ ■ , Qf^- d — min(n + 1, r + 1) since Q'f^ = if 
a>r + 2. 
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coefficients {bn\o)} up to the order of 5, where we put A = cosh?]. 

b?iO) = '-^, (4.2) 

6A^Qf' (2A^ + l)Qf^ (4A^-l)gf 
+ ^a)2 + 



4.(0) . -^J^ + l™^ + ilf^^, (4.3) 

4..(o)^gggg^,ggi!|!,(^ + ^t^'^'' 

, (8A - 32A') Q'^'qf ' , (-12A^-6A)Qf'' + (8A'-4A)QS-' 

Qi'^ Q? 

(1) 140A^gf^' , (240A^-60A^)gf^gf^' 

"4 {^) - 7778 ^ 7777 

(^_14A4 - _ 1^ g(2)' + (24A4 - SA^ - 1) gj^) 

(80A^ + 40A^) gf + (40A^ - 80A^) Q^^'^Q? 
(-40A4 + 20A2 - I) gf 

(_96A^ - 8A^ + 4) gf^gf ^ + (16A^ - 12 A^ + 1) gf 



(4.5) 
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5Q\" gr 

+ ¥ + f ) ^ (8A^ + 10A3 + f ) gf ) 



2 



-12A5 - 30A3 - 8A) gf^ + (40A5 - 6A) gj^^ 



(1)4 
1 



-560A5 - 280A3) gf + (672A5 - SSGA^) Q^^Qf^ 



(672A5 - 336A3 + 42A) gS^^^gf ^ 



(-160A5 - 100A3 + llA) gf ^gf ^ + (64A5 - 40A3) gf ^ 
(960A5 + 120A3 - 60A) Qf^Qf^^ + (-192A5 + 144A3 - 12A) gf ^gf ^ 
(-192A5 + 144A3 - 24A) gf ^gf ^ 



fm^ + + 20A) gf + (-320A5 + 80A3 + 30A) gf ^gS' 



(4)^(2) 



(40A3 - 160A5) Q\'> + (32A5 - 32A3 + 6A) g^' 



^ — — . (4.6) 

.(1)6 ^ ' 



In deriving these coefficients ()4.2|) - ()4.6|) . we did not assume ()2.40|) . Of course, 
when one calculate the free energy of the model, one must assume ()2.40j) and 
()2.48p . We can also rewrite the coefficient bn^ (0) in terms of gf gf • • • , gf ■* 
and ( { d = min(?T, + 1, r + s + 1) ) since gf for a G Z>r+s+2 can be written 
in terms of Q^i\Qf\ ■ ■ ■ , gf and ( due to the relation ()2.3()|1 in the 
limit V irj^^oo (see also an example: ()2.4H) - ()2.43j) ). If 6f''(0) is written in 
terms of gf Qi'\ ■ • • , gf and C, { d = min(n + 1, r + s + 1)), it should be 
the coefficient of the high temperature expansion directly derived from the 
reduced NLIE ()3.9|) - ()3.17|) . Of course these two expressions of the coefficient 

^Qf = C if s = -1. 
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0.5 1 1.5 2 2.5 3 

T 



Figure 1: Temperature dependence of the high temperature expansion of the 
specific heat C for the rank 2 case (r + s = 1, J = 1, g = 1, = (a G B)). 
We have plotted plan series (dotted lines) of C in Appendix and their Pade 
approximations of order [n,d] (numerator: a degree n polynomial of 1/T, 
denominator: a degree d polynomial of 1/T) by using Mathematica: each 
line denotes C for s/(3|0) with [20,20] (thin), s/(2|l) with [17,17] (medium), 
s/(l|2) with [17,17] (thick), s/(0|3) [20,20] (dashed thick) respectively. We 
have also plotted (thick dots) a result of numerical calculation from another 
NLIE by Jiittner and Kliimper [52] for the s/(2|l) case. C for the s/(3|0) case 
was also considered in |57l 16] . 

6n''(0) are equivalent under the relations ()2.40p and ()2.48p . 

For fixed values of parameters, we have calculated the high tempera- 
ture expansion for much higher order (see, appendix). We have plotted the 
high temperature expansion of the specific heat (Figure ^IHI). Here we have 
adopted the Pade approximation method. There is a duality among the 
specific heats with respect to interchange of r and s. In particular, r = s 
case is self-dual, then the specific heat becomes an even function of T (see 
()A.1.6j) ). In Figure CJ we have also plotted a result of a numerical calculation 
by another NLIE [53j. We find a good agreement between our result and 
their result except for very low temperature region. 

We can also calculate the high temperature expansion from the NLIE for 
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Figure 2: Temperature dependence of the high temperature expansion of 
the specific heat C for the rank 3 case (r + s = 2, J = 1, g = 1, = 
(a G B)). We have plotted plan series (dotted lines) of C in Appendix and 
their Fade approximations of order [n,d] (numerator: a degree n polynomial 
of 1/T, denominator: a degree d polynomial of 1/T): each line denotes C for 
s/(4|0) with [19,20] (thin), s/(3|l) with [17,17] (medium), s/(2|2) with [16,16] 
(thick), s/(l|3) with [17,17] (dashed medium), s/(0|4) with [18,21] (dashed 
thick) respectively. C for the sZ(4|0) case was also considered in 
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Figure 3: Temperature dependence of the high temperature expansion of 
the specific heat C for the rank 4 case (r + s = 3, J = 1, g = 1, /!„ = 
(a e B)). We have plotted plan series (dotted lines) of C in Appendix and 
their Fade approximations of order [n,d] (numerator: a degree n polynomial 
of 1/T, denominator: a degree d polynomial of l/T): each line denotes 
C for s/(5|0) with [17,21] (thin), s/(4|l) with [16,18] (medium), s/(3|2) with 
[17,17] (thick), s/(2|3) with [16,17] (dashed thin), s/(l|4) with [16,18] (dashed 
medium), s/(0|5) with [17,21] (dashed thick) respectively. 
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= — 1 in subsection 3.2. Similar to ^ = 1 case, we assume 

(^''3^ J J \ 

r«(t;) =exp Y.^mAv){7^T + 0{{^Y^'^') , (4.7) 

\n=0 / 

where hmfl{v) = logQm''. Here we do not expand {Qk^^}k>i with respect to 
^. ()4.H1 for a = 1 should coincide with ()4.7j) for m = 1 up to a factor from 
the normalization function ()2.17|) . Thus we have 

&i'H0)=6i,„(0) + 2A5„,i (4.8) 

Due to symmetry between the NLIE for ^ = 1 and the one for ^ = — 1, the 
following relation follows: 

6,„(0) = (-l)"6«(0)|^,„,^^,),^^^^,. (4.9) 

For example, ()4.8|1 and ()4.9j) for n = 1 and ()4.2p reproduce the Q-system 
(jm)) for (a, m) = (1, 1). From the relations (jOj) and (ji^ for n = 2 and 
()4.3|) . we obtain identities among characters 

-3Qf + )g«' + 2gf )g« = -sg^' + g«g«' + 2g«g«, (4.io) 

gf^gS^^-gf^ = g^^^gS^^-g^^\ (4.11) 

where we have used the fact that Q-m does not depend on A. These relations 
can be proved from the relations ()2.25p . ()2.2(j|l and ()2.4()p . 

Some comments on references on the high temperature expansion are in 
order. The high temperature expansion of the free energy was calculated 
from the Takahashi's NLIE for the XXX-model up to the order of 100 
|46j : the XXZ-model up to the order of 99 JH]- As for the higher rank 
or higher spin case, we have some results fHl ^1 ^1 HB] from NLIE. In 
particular, our result on the s/(r + 1) Uimin- Sutherland model in was 
applied jl?! EHl EHl EHl HH Ii2| to spin ladder models and good agreement was 
seen between theoretical results and experimental data. We note that the 
coefficients ()4.2j) - ()4.4|) coincide with eqs. (4.14)-(4.16) in ^Hl- Note however 
that the coefficients in our paper are more general than the ones in ^H] since 
the value of Q[^^ flCTIl was restricted to s = — 1 case in ^Hj- There are also 
several works on high temperature expansions by different methods (see for 
example, |S1 ESI ESI EZl EH!)- 

5 Concluding remarks 

In this paper, we have derived NLIE which contain only r + s + 1 unknown 
functions for the U{sl{r + l|s + 1)) Perk-Schultz model. The key is a duality 



24 



for the auxiliary function ()2.23p and the quantum (supersymmetric) Jacobi- 
Trudi and Giambelh formula ()2.25|) and ()2.26|) . Although we assumed that 
q is generic, we expect that our NLIE (at least reduced ones ()3.9|) - ()3.17|) . 
()3.2H) - ()3.29|) ) will also be valid even for the case where q is root of unity 
as we will not need to take into account truncation of the T-system. The 
high temperature expansion of the free energy in terms of characters was 
calculated from our NLIE. 

There are NLIE with a finite number of unknown functions for algebras 
of arbitrary rank in different context [SHI lEIl • These NLIE are different from 
Takahashi-type. Whether one can generalize (or modify) their NLIE for finite 
temperature case is still not clear. A deeper understanding of this subject is 
desirable. 

There is an another kind of formulation of transfer matrices which is based 
on the graded formulation of the quantum inverse scattering method. In this 
formulation, the row-to-row transfer matrix is defined as a supertrace: t{v) = 
stro{RoL{v) ■ ■ ■ Ro2{v)Roi{v)), where the i?-matrix is defined as Ral'lli"^) ~ 
(^—l^pCaOpCf-O/^^i'^i^^) and the graded tensor product is adopted. As far as 
the free energy (in the thermodynamic limit) is concerned, we think that 
there is no difference between this graded formulation and the one we have 
adopted. 
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Appendix: The high temperature expansion of the specific heat 

We will list the high temperature expansion of the specific heat Csi{r+i\s+i) for the Uq{sl{r + l\s + 1)) Perk-Schultz 
model at g = 1, /ia = (a G B). Here we put t = In this case, Q^^^ (cf. ()2.40|) ) becomes 

j=0 



which is the dimension of a-th anti- (super) symmetric tensor representation of sl{r + l|s + 1). If one substitute 
flA.l.l|) . A = 1 and the values of (r, s) into fl4.2|) - ()4.6p . one can recover ()A.1.2|) - ()A.1.9|) up to the order of 5 through 

C = —Tq^. a formula for r < s can be obtained from the relation Csi{s+i\r+i) = Csi(r+i\s+i)\t^-t- 

8t^ 16t^ 40t*' 400t^ 1246t'' 11228t'' 43343t* 649298*^ 120769*^" 5559367*^^ 36953579*^^ 4333458857*^^ 298222277909*^* 

to '^sUalo) ^ 9 + 27 27 243 729 ^ 3645 32805 137781 918540 ^ 885735 ^ 18600435 584585100 58926178080 

44130279500393*1^ 2885993845291237*^'' 47755357995530701*^'' 4655618035381741733*^^ 45436230289647581*^'' 1590783575674338541*^" 

+ + + + 

5745302362800 321736932316800 7239080977128000 347475886902144000 12306437661117600 89350942346265600 

11365203602766081451*^1 1297481476315241042696509*^^ 8454419484929269090011049*^^ 46816780786984484371594673*^*1 8261193033376436054715299*^'' 

H + + 

8102505908218176000 60606744193471956480000 954556221047183314560000 2016022738851651160350720 445851182630653622000640 

4614044757865223484648570791543*^'' 16687567624201209926686045552339*^''' 764945949570840761119882192959107*^'* 4370526181353809390696401823321023*^" 

h 1 

208658353471145895096299520000 558906303940569361865088000000 45209309918748277270864896000000 104627260097674584541144473600000 

4077089856720735402715997482797183733*^° 469645986012529902658517386886308977221*31 5501623615004327359193974711230889260281*^^ 

615208289374326557101929504768000000 8920520195927735077977977819136000000 583888594642542659649467639070720000000 

81517987350487140844545467182506851908591*'''' 8539638490748569692670776190970340550336847*3'1 565870839129464697660769748292242672448332479*^^ 

+ + 

1351398263165884933985080078663680000000 273059671917403950089785323323129856000000 9097613107632737375587559089563893760000000 

2803571976313303389947366028586799714153992385183*3'' 821880309698434533395036032147535806012394330483*3'' 

1 

48253739922884039040116413411046892503040000000 14814744713166152336877846222689835417600000000 

58773945021047530582522114436884890912882464905123*3* 903065874685632945085489557823621 181308117028891323*3" 

+ 

668128706624548232863150339537572357734400000000 24102743091480577500538148498817922805268480000000 

69053384918361487529760006169534549420582996627140401**"' 

+0(**1) A. 1.2) 

586849397009961886969624485188610294389145600000000 



32t^ 3()4t-' 5480t* 8320*'' 736708*" 1470644t^ 146834891*** 14149151840*'' 228231260059*^" 3300909899909*11 86046211353427*1 

+ . ■ ■ ■ ■ 



■' ' 27 729 2187 6561 177147 531441 23914845 2711943423 27119434230 366112362105 7908027021468 

6841741391685967*1^ 41669257450473325*1* 209509171518293955313*1^ 23504662660033768183787*1" 59568189209735825524303*1'' 



466008736193660 3131578700501328 9159867698966384400 1B388B77734263B2679200 1731214996104646651600 

361017420632075630992718067*1^ 3424600460806080368078749*1^ 8630572663979741453598354479*^" 17202731244586324123474774048501*^ 



22436546336556220604736000 68110944235974241121520 588478558198817443289932800 240139375283176527142516896000 

74148256847472328976368306828013*^^ 10062652817270791187839874880710933*^^ 708322746602409944512187717871912169*^* 



7924599384344825395703057568000 100858537618934141399857096320000 316324648705015526355199808330342400 

15024962981052794085153555040621422973*^'' 231254246428418304396673721279263409933713*^" 38713134042756953404786472697685030409352013*^''' 



110457493727965947968161876390656000 9821880342290732093328954048657131520000 213099725283629276667762128019971692800000 

174322192700322839649823295381972425935690253*^** 15926771354416351317222940636770636517312076257*^" 



2938193181940949117691871765123852128000000 67019011202800272994904518213769017498828800 
668493937475506754534088912524547499759600367865761*^° 336260111063276140660586885786196146564962534202441623*^ 



5757542326058750725471342701091974685126656000000 1101993601207644888865214992989003964733241968400000 
24356652704738251987509066739814466042906238637219292447*^^ 238144042414431030672413060665368221701372461889917292229*^ 



120217483768106715147841635598800431425444577280000000 626042066206424505863222816152802925415196958720000000 



19283712635849148986736121736756878449365536043407449116181f 
58382894402793416775358836340649964244434367807488000000 



34 



+ 0(t^'') (A.1.3) 



15*^ 15*^^ 435** 345*'' 9555*" 21917*'' 172967*** 3052445*^ 53684587*1" 41381153*" 24190901579*1^ 74629743461*1^ 59210969497* 
Csi(4|0) ~ + 32 256 256 4096 8192 65536 688128 22020096 6291456 



6291456 15854469120 8304721920 186025771008 

831873403828903*1^ 40380622501051099*1" 69674366936531941*1'' 677763076075244557*1* 27733137112330033541*1" 916368739307996439457*^" 



72550050693120 12188408516444160 5078503548518400 88642971028684800 1808316608985169920 68199369253154979840 

195301776246305171789377*^1 71292804015129538500833*^^ 2200476451580705384154142457*^^ 78446060458744023170123124563*^* 



12368885606458562252800 3445590165496528896000 152384670659249486954496000 2681970203602790970399129600 

125238936013261236946687178727*^^ 191117971515476319277466700295697*^" 1065649288731276091956499644404357*^'^ 



11862560515935421599842304000 4934825174629135385534398464000 317238761226158703355782758400000 

11192351020362615444556022460039019153*^* 156907755522484105311278338009333117*^" 2360085804850806890473467996881990318573*^ 



230949818172643536043009848115200000 19795698700512303089400844124160000 41081897067886708999650692982374400000 

29423655375616416313737080966732640595067*-'*1 3157300240979374909778816498331 173586672649***^ 16704208032661485316059904933586012489234849*** 



1227477288149584699201684341837004800000 49099091525983387968067373673480192000000 369326254640301513906859729911545856000000 

4632447279343820265500763653186467289421668361*** 61387162556304338047936313287382785565501241**^ 



688846225797682366509708661 9B()734B7()848()()()0() 8550750332321 1849270925129324009881 OOOOOOO 

!n74572K524i44()73:!8571()4S:jS48y(j:!7:j:i2:!lj2:i8i)241lj:!Srji:!/'"'" i:!7i5:!5i)i!)i2i5924(i847()8:!7rj 1 i!)uy()!)U2y 72:!()(i7S<J()48i 78:! 



21640793029659989226269007324158440419360768000000 13288206246282449524902022041149919555747840000000 

20934050741656543039384025455100076026671071927370779*** 10632927999476868411724450072247106687260120276680120833**" An 

+0**" A.1.4 

399622332856261339661889365984463615434352640000000 76868096841352281731707514015410780369669447936000000 



(i9t^ 45t^ 8745t* 4()65f'' 4474()5t'' 2734683*'' 406592713*** 222295155t^ 359912058803(1" 2122554602333tll 143332776011113* 
a!(3|l) = — — + 



128 4096 4096 131072 1310720 83886080 58720256 56371445760 338228674560 18038862643200 

2496231619276031*1^ 2384568827915515*1'' 4333850790231468637*1'* 13960280493348579178073*1'^ 6984410297833152298633*1'' 633205038656776496328727*1 



266121057382400 263987186016266 297165007639019520 1331299234222807449600 332824808656701862400 58093057493358870528000 

14162457941778109685145689*1^ 604540593802926375008290949*^° 1296710847223129254977218906963*^1 86141375198512658190978194262097*^^ 



482817866611471501721600 59593518178330196783923200 32424291481673293431980032000 11413350601513799288056971264000 

99308058853685747403471990563649*23 49493071668124032061338672924297499*2" 6931634160144673731055617563665688207*25 



1862318279967286597118853120000 22497996706704001166617901766596800 99510370044460006115809950072832000 

587874112017235875229135097766679692569*2'^ 29312181553180719038186824685428443001969*2'' 1698158664136245991126625182277922005467823563*2 



82792627876990724256353878460596224000 328542174115042556572832851034112000000 77493899692863317903947230239118065664000000 

3126686061624145346779158233961466863595139419*29 197403265536059174349045556856630021916456067183*3° 



27897803889430794445421002886082503639040000 4463648622308927111267360461773200582246400000 
3751792081077294678466899166499824476489567047662019*31 80680001031677274734626646657250828568042366729432743*3 



27183620109861366107618225212198791546880676000000 1054394961837046927810646311261044035719004160000000 
1319799435656923687276493800076182418349748116331178161*33 6324243289350644149071369870231560459419389643237919619*'' 
7931220926171316229046258660147412121621626880000000 51904772136387320644826041672092637674131308644000000 



(A. 1.5) 



9*2 305**1 4165*^ 1028409*** 758654369*1" 775100578187*12 2108183654669*1* 36086372927030207*1'' 123454173470857039087*1* 



^ '=s!(2|2) ^ + 163840 ^ 82575360 59454259200 ^ 116266106880 1451001013862400 ^ 3656522554933248000 

775360975454089529227*2'^ 414116620493362593763666669*22 18120085561666997913793728601*2* 4422094856669703488323197127729t2'3 



17049842313288744960 6802887083002209239040000 223497516966899247533260800 41123543121909461546119987200 

61544192334286763277931264839087079063*28 24624774786846908325218781376665239554743*30 1066449426872776917417866273408875411974943*3 



433030909073706630080643465216000000 130948546903888884936386583881318400000 4332272781704416585417709441777664000000 

3971063172619299668637043400605041165950069531*3'' 

+0(t3S) (A.1.6) 

12300825460672899401959083249298833408000000 



24*^ 48t3 1128t* ;-i504f'' 8274*'' 178444*'' 1306457*** 160692418*'' 3091451869*1" 177674519*^1 173473316029*1^ 4227122268577*1 

+ 



25 125 625 3125 3125 78125 390625 41015625 820312500 29296875 46142578125 483398437500 

421192483420837*1* 157545290437200577*1^ 499112291761171129*1'^ 20544397967491432423*1'^ 2622947397527190428811*1* 



136361662600000 13196777343750000 316722656250000000 1319677734376000000 2217068693760000000000 

344106766804640679764433*1^ 6481099627451266704610167*^° 514322825251322335204973701*^1 8742897674113739154665330649*^ 



17667185668945312500000 1211464160166260000000000 21971664641016626000000000 767260634765625000000000000 

11410614110336974477848290940603*^^ 3647621524555860483886516418494301*^* 6101990866174335716593365406813719*^^ 



422962424914660781250000000000 186099066962402343750000000000000 206782622121887207031260000000000 

517943244817437855867592189832638271*^" 21134646896102733105412988229490963037*^'' 9912970564455830832161043262688106009427*^* 

+ + 

17121114160541015625000000000000000 687901908236022949218750000000000000 227632995089011230468750000000000000000 

1100169921280159675483593435467841886949801*^" 14414030290129165020160816003617622620734617*3" 11096391762725899024492170038454077512197326053*^1 



37569444189686853027343750000000000000000 242071098102329589843750000000000000000000 467474030230385869873046875000000000000000000 

1298913266483771488098839606916528516217512843843*^^ 36827816163436004340581866161243966139756265512563*^^ 



16636419281104940722666250000000000000000000000 2502652047730934464599609375000000000000000000000 
5904769218118219660609568211766656164631746178169*^* 93947676821931398956438041874300826546727311316769613*^ 



69827681730247639062500000000000000000000000000 46799593292568474488012696312600000000000000000000000 
6904178434667906871447348617674173162386567620702330884299*^^ 100666404268726322722978537972143488460344388408692120437*^ 



67282702190103812773327639062500000000000000000000000000 3823216202619438304866957031250000000000000000000000000 
286059137038561139282916743342892323996039400157684633307181*^* 



to 2011699909686919846654980468750000000000000000000000000000 

CO 



- 0(f"') (A.1.7) 



CsiC-lll) - 



+ 



912t'' 1591764* 8()87472f°' 1583091021" 10619486948*'' 139915063629** 50193776068378t'' 145147029308647729*^" 



625 3125 78125 9765625 48828125 6103515625 30517578125 16021728515625 24032592773437500 

154561110991941973*^^ 511262973476333776477*^^ 562798123639186677983*^^ 27263166260140540122161*^* 72774085086186213945892331*^ 



30040740966796875 67591667175292968750 70810317993164062600 3004074096679687500000 619B9028244018BB4687B0000 

26166239133667469488082468441*1'' 2149440171168061134748079150611*1'' 11616816820943999916517438797466393*1* 44866031236953872796448471717341173*1^ 



2606898476646972656250000000 128142B3B68649291992187B000000 1014888882637023926781260000000000 192BB70424646139144897460937B00000 

272821227014028414160064869868461716981*^° 66246387219780884484696724793606687399357*^1 266119119131347681110777846066146382444833999*^^ 

+ 

23106846096763669738769531260000000000 2096370726728671414947509765625000000000 23049077983014286664422607421876000000000000 

8476891798073835488652808808090088670094654917*23 434354342487729446922438313016331279852618182267*2'1 



201679432351374998688697814941406250000000000 49299416797002777457237243652343750000000000000 
74214528639575465586270534190435386898519833666641*2^ 241489502492350466125582435489379440678333756529927*2" 



1352823539616565540200099349021911621093750000000 57851356445462442934513092041015626000000000000000 
360758243689593696611421633292891289316448278820393793389*2'^ 209490702367563791064096969866392414710044648636093273111309*2* 



6126268610090188303729519248008728027343760000000000000 55967933222184856276726350188256310058593750000000000000000 
372970956002091736988968941025697451968680727107686289984993771*23 4270970709112276963496423117872370698202654386252499239578642356083** 



4197594991663864220764476264119148264394531260000000000000000 264448484474823446907532004639506340026866468760000000000000000000 
705454893266763100766781948219461396410990670736178037506736499754037**1 40167992980468192651166832646604563149328974033306535140279044661346513**2 



6390838374808233276098690112121403217315673828125000000000000000000 1161970613601496959290670929476618766784667968750000000000000000000000 
137948929879317690635779621721238444356037066949039424111273171238709303*** 



CO 1022270328271698821254176436923444271087646484375000000000000000000000 

295393994896927994224867184092844725554 7963144680752974314176599201612597751*-^* 

hOr**^) (A.1.8) 

48518201539360464871606382075697183609008789062600000000000000000000000000 



672*^ 336t^ 1722964* 3045504*^ 1792271881" 4080719804*'' lfi814157()529t** 198377fi738r)21fit'' 1492774189466571*1" 63203192343466109*1^ 
'^sl(3\2) = 1 1 1 ^ H ■ 



625 3125 78125 9765625 48828125 6103515625 30517578125 16021728515625 190734863281250 30040740966796875 

484746376821625940723*1^ 119657043279360047669*1^ 34147743686560410475681*1"' 20015701333284035643585049*1^ 1296881833233097980209112581867*1" 



45061111460196312500 36406168996682031260 2360343933106468760000 383BBB8891296386718760000 67659258842468261718760000000 

306368982994386261087366331797*1'^ 1412566068743727245477698767563469*1* 38642600881761570643294936210984427*13 



39164663681983947763906260000 56382715702056884766625000000000 3369748243130743503570666640626000 

749943826781288708352366220297611952081*^° 823393490024196742631909258671419921007*^1 30040873842462516099488089062008818909101223*^ 



23106846096763669738769631260000000000 49889779184013605117797861562600000000 720283686969196423888206481933693760000000 

47197534334910699474197303571 154961451 19402591*^-'* 1816573282828621256378963795157922371 189821325851*^"' 



2016794323513749986886978149414()625()()00()0()00() 34130365474848076701164245605468750000000000000 
8322591 197319584873907765768566656873942892240011*^'°' 78284346134584453501416727068291711888925901387282239*^ 



253946171687857713660329208374023437500000000000 1159656736020406242460012435913085937500000000000000 
58189014600621207023093446380939189933968282293477470193*^'^ 282847733617417987601884583482566260456466391798062174307*^ 



1281317162522647075932379812002182006836937500000000000 33219333B839178871B38024396704269409179687B000000000000 
130877066690281430610636720696462523358644877286866176351960231*^3 4709049750740948620371923839229876772436358649661968355418846451943*^ 



209879749B831932110377238132069B7412719726562B000000000000000 44074747412470674317922000773251066671142678126000000000000000000 
11079182595131990060363516685430527467174991034289759600690770327879*^1 51676596103109262978856677266778347122077660630622361028974962093289069*^ 



130426272966270066869156941063702106476830078126000000000000000000 387323637867166653096890309825539588928222656250000000000000000000000 
1673903716406893043636199894308218634748017092628006510068625803350703093*^^ 



CO 1456735217787170820287201422(il59()8()S62998962402343750000000000000()0000() 

2560112840573273875628881469331020252804964139812191437239808660979503389001*3* <,,; 

+ +0(t3^) (A.1.9) 

16437609680706602469147485206903649330139160156250000000000000000000000000 
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